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Abstract 

In the paper stochastic Volterra equations of nonscalar type in Hilbert space are studied. The 
aim of the paper is to provide some results on stochastic convolution and mild solutions to those 
Volterra equations. The motivation of the paper comes from a model of aging viscoelastic materi- 
als. The pseudo-resolvent approach is used. 
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1 Introduction 

Assume that H is a separable Hilbert space with a scalar product {-,-)h and a norm \-\h- 
We are concerned with a stochastic linear Volterra equation of nonscalar type in H of the 
form ^ ^ 

X(t)=X + [ A(t-r)X(r)dr + [ ^{r)dW{r), t € [0,T\. (1) 
Jo Jo 

Let G denote another separable Hilbert space such that G =— > H (G is densely imbedded 
in H). In (JTJ), A € L 1 1 oc (R + ; B(G,H)), where B(G,H) denotes the space of all bounded 
linear operators acting from G into H. We assume that W is a cylindrical Wiener process 
with values in some another Hilbert space U, and ^ is a stochastic process specified below. 

Observe that Volterra equation of so-called scalar type (see (Priiss 1993)) corresponds 
to the operator A(t) of the form A(t) = a(t)A, where a € (R+), and A is a closed linear 
unbounded operator densely denned in H and G = (D(A), \-\a), where D(A) denotes the 
domain of A equipped with the graph norm | • \a of A. 

The setting for Q when A(t) are bounded seems to cover many abstract treatments 
of Volterra equations although this is not the most general approach. However, it is very 
difficult to develop a reasonable theory in very general case even for deterministic equation. 

As we have already written, the paper has been inspired by a model of a hereditarily- 
elastic anisotropic aging body with a straight crack, see (Costabel et al. 2004). 

The equation Q is a stochastic version of the deterministic Volterra equation 

u(t)= [ A(t-T)u(r)dT + f{t) , t>0, (2) 
Jo 

studied extensively in (Priiss 1993). 
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Definition 1 A family (S(t))t>o C B(H) is called pseudo-resolvent for J1J) if the fol- 
lowing conditions are satisfied: 

(51) S(t) is strongly continuous in H on K + , and S(0) = I; 

(52) U(t) = J *5(r)dr is leaving G invariant, and (U(t)) t >o C B(G) is locally Lipschitz 
on R + ; 

(53) the following resolvent equations hold 

S(t)y = y+ [ A(t-T)dU(r)y, t > 0, y G G , (3) 
J o 

S(t)y = y+ f S(t-r)A(r)ydr, t > 0, y G G . (4) 

JO 

• Equations 0) and (0) are called the first, resp. second resolvent equation 
/or H. 

• ^4 pseudo-resolvent S(t) is called resolvent for if * n addidion 

(54) for y G G, 5'(-)y G G a.e. and S(-)y is Bochner-measurable in G on R_|_. 

Comment Because A(t) G B(G,H) the convolution A*dU in |3p is well-defined though 
the function A{-)y is not assumed to be continuous. Indeed, for y G G i/ie function 
g{t) := U(t)y is locally Lipschitz in G by the condition (S2) and then g G BV\ oc (M + ; G). 
For details, see (Priiss 1993, Section 6). 

Let us emphasize that the pseudo-resolvents are always unique. We shall assume in the 
paper that the equation Q is well-posed in this sense that Q admits a pseudo-resolvent 
Sit), t > 0. Precise definition of well-posedness is given in (Priiss 1993). That definition 
is a direct extension of well-posedness of Cauchy problems and Volterra equations of scalar 
type. The lack of well-posedness leads to distribution resolvents, see e.g. (DaPrato and 
Iannelli 1984). 

In sections El and 01 we shall use the following Volterra Assumptions (abbr. (VA)): 

1. AeLl c iR + ;B(G,H)); 

2. S(t), t > 0, are pseudo-resolvent operators corresponding to (^). 

2 Probabilistic background 

We are given a probability space (fi, T, (^i), P), t > 0, with normal filtration and a 
Wiener process with the positive linear covariance operator Q in U. Assume that W 
is a cylindrical process, that is TrQ = +oo (if TrG; < +oo then W is a genuine Wiener 
process). We will need the subspace Uo := Q l l 2 {U) of the space U, which endowed with 
the inner product (u,v)u '■= (Q~ 1 / 2 u,Q~ 1 / 2 v)u forms a Hilbert space. 

A linear, bounded operator B acting from Uq into H is called a Hilbert- Schmidt if 
Ylt^i \Buk\ 2 H < +oo, where {u^} C Uq is a base in [To- The set L2(Uq,H) of all Hilbert- 
Schmidt operators from Uq into i7, equipped with the norm \B\£ a m jj\ := (X^fcS l-^fclij) 1 ^ 2 ) 
is a separable Hilbert space. For abbreviation we shall denote L\ '■= L2(Uq, H). 

Assume that ^ belongs to the class of L^-predictable processes satisfying condition 

p(^ T |^)lioaV<+oc) =1. 
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Such processes are called stochastically integrable on [0, T] and create a linear space which 
will be denoted by M(0, T; L° 2 ). 

Let &(t), t G [0,T], be a measurable Lg-valued process. Define the norms 



||$||t : = jEQVMlio^)} 2 = I® J* [lH*(T)Q3)(*(r)Q3)*] drV , t€[0,T\. 

By M 2 (0, T; L®) we shall denote a Hilbert space of all Lg-predictable processes 3> such 
that ||$||t < +oo- 

In the whole paper we shall use the following Probability Assumptions (abbr. (PA)): 

1. Xq is an H -valued, J-q -measurable random variable; 

2. <f G A/" 2 (0, T; Lg), w/iere i/ie interval [0,T] is fixed. 

Definition 2 Assume that conditions (VA) and (PA) hold. An H -valued predictable pro- 
cess X(t), t G [0,T], is said to be a strong solution to (QJ), if X has a version such that 
P(X(t) G G) = 1 for almost all t G [0,T], J Q T \A(t - T)X(T)\ H dT < +oo, P - a.s. and 
for any t G [0, T] the equation (QJ) holds P-a.s. 

Definition 3 Let conditions (VA) and (PA) hold. An H -valued predictable process X(t), t G 
[0, T\, is said to be a mild solution to the stochastic Volterra equation 0), if 

E(J\S(t-T)^(T)\ 2 L odr]<+oo for t<T (5) 
and, for arbitrary t G [0, T], 

X(t) = S(t)X + f S(t - r)*(r) dW(r), (6) 
Jo 

where S(t) is the pseudo-resolvent for the equation (0). 

Comment The well-posedness of £7J) implies existence and uniqueness of the resolvent 
S(t), t>0, and in the consequence, existence and uniqueness of the mild solution to the 
equation Q^- 

Proposition 1 Assume that (VA) hold, B is a linear bounded operator acting from U 
into the space H and 

[ \S(r)B\ 2 L odT= [ Tt[S(t)BQB*S*(t)] dr < +oo, (7) 
Jo 2 Jo 

where B*,S*(t) are appropriate adjoint operators. 
Then 

(i) the process W B := £ S(t-T)BdW(r) is Gaussian, mean-square continuous on [0, T] 
and then has a predictable version; 

(") 

CovW B (t)=[[S{T)BQB*S*{T)]dT, te[0,T}; (8) 
Jo 
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(iii) trajectories of the process W B are P-a.s. square integrable on [0,T]. 

Comment Proposition^ is analogous to the known result obtained in the semigroup case. 
Gaussianity of the process W B follows from the definition and properties of stochastic 
integral. Part (ii) comes from theory of stochastic integral while (iii) follows from the 
definition of W B and the assumption 

Proposition 2 Assume that conditions (VA) hold with operator A G W 1 ' ([0, T];B(G, H)). 
Let X be a strong solution to the equation (QJ) in the case *&(t) = B, where B G B(U, H) and 
trajectories of X are integrable w.p. 1 on [0,T]. Then, for any function £ G C 1 ([0, T\; H*) 
and t G [0, T], the following formula holds 

(X(t),£(t))H = (X o ,C(0)) H + f\(A*X)(T)+A(0)X(TU(T)) H dr 

Jo 

+ f (^r),BdW(r)) H + ! (X(r),i(r)) H dT, (9) 
Jo Jo 

where dots above A and £ mean time derivatives and * means the convolution. 

Proof First, we consider functions £(t) := ^o^C 7 ") ^ r £ 1 where £o S H* and 

<p G C X %T\. Denote F^{t) := (X(t),£ ) H , t G [0,T\. 
Using Ito's formula to the process F^ (t)(p(t), we obtain 

d[F^(t) V {t)] = <p(t)dFz (t) + tp(t)Fz (t)dt, t G [0,T]. (10) 

Then 

= {f A(t-T)X(T)dT + A(Q)X(t),£ ) H dt+(BdW(t),t ) H 
Jo 

= ((A*X)(t)+A(0)X(t),€ o ) H dt+{BdW{t),Z o )H. (11) 
From and (fTTj) . 



F €o (t) V (t) = F&(0M0)+ V (a)((A*X)( a )+A(0)X(a),e >ff<fa 

JO 

+ / M«)Bd^(s),^ >H + / <p(s)(X(s),Zo)Hds 
Jo Jo 

= (X o ,£(0))+ [ ({A*X)( 8 ) + A(p)X(8),t(8)) H ds 
Jo 

+ / (BdW(s),tts)) H + [ (X( S ),i(s)) H ds. 
Jo Jo 

Hence, we proved the formula @ for functions £ of the form £(s) = ^oV 9 ! 5 )) s £ 
[0, T]. Because such functions form a dense subspace in C 1 ([0, T]; -fP), the proposition 
is true. ■ 
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3 Properties of stochastic convolution 

In this section we study mild solution to the equation (^Q). We use the following notation 

W*(t):= fs{t-T)nr)dW{r\ t€[0,T\, (12) 
J o 

where the condition Q holds. 

Proposition 3 Assume thatS(t), t>0, are the pseudo-resolvent operators to the Volterra 
equation P| ). Then, for arbitrary ^ G Af(0,T; L®), the process W^(t), t > 0, given by 
has a predictable version. 

Proof Proof is analogous to construction of stochastic integral, see e.g. (Liptser and 
Shiryayev 1973). 

The process S(t-r)^(r), r G [0,T], belongs to J\f(0, T; L%), because * G J\f(0,T;L%). 
We may use the following estimate: for arbitrary a > 0, b > and t G [0, T], 

PQW*{t)\ H > a) < 1 +P \S(t - r)^(r)\lodr > ft) . (13) 

Since pseudo-resolvent operators S(t), t > 0, are uniformly bounded on compact itervals, 
there exists a constant C > 0, such that \S(t - r)*(r)|| < C 2 |^ / (r)|^ , r G [0,T]. 

Then (JEU) reads 



P(|W*(t)|H >a)<^ + p(J* \Hr)\lodr > ± 



(14) 



We prove predictability of the process in two steps. In the first step ^ is an 
elementary process, so the process W has a predictable version by Proposition^ part (i). 

In the second step 'F G J\f(0,T; L^). There exists a sequence (^ n ) of elementary 
processes that for arbitrary c > 0, 



P (£ Mr) - Mr)\ 2 L odr > c) ^ . 



(15) 



Because, by the previous part of the proof, the sequence (W*) converges in probabil- 
ity, it has subsequence converging almost surely. This fact implies the predictability of 
t G [0,T]. ■ 

Proposition 4 Assume that * G AA 2 (0,T;L^). T/ten ffce process W*(t), t > 0, denned 
6y fjjjj) /ias square integrable trajectories. 

Proof From Fubini's theorem and boundness of operators S(t) we obtain 



E 



o 



S(t - t)^{t)oIW{t) 

7 

JO 



|5(i -r)*(r) | 2 drdi < M 



E 
T rt 



S(t - T)^{T)dW{T) 







H 



dt 



o Jo 



\V (t)\ 2 l o drdt < +oc, t€[0,T). 
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Theorem 1 Assume that (VA) and (PA) hold with A G W 1 ' 1 ([0,T],B(G, H)), and X is 
an H -valued predictable process. Let S(t) be a weak analytic resolvent for 0). Then strong 
solution to (QJ) is a mild solution, that is 

X(t)=X + W i!! (t) tE[0,T]. (16) 

Proof By Proposition^ for any £ G C^QO, T], H*) and i G [0, T], the following equation 

holds 

(X(tU(t)) B = (X ,Z(0))H+ [\(A*X)(T)+A(0)X(T),Z(T)) H dT (17) 



+ / (^(T)dW(T),C(T)) H + / (X(r),£(r))^dr, P - a.s. 



Let us take f(r) := S*(t - r)C, for ( G t 6 [0,t]. 
Now, the equation (j!7j) may be written like 

(X(t),5*(0)C)i/ = (X ,5*(t)C)i/+ A(i*X)(r)+A(0)X(r),5*(t-r)C)ffdr 



+ / (*(T)dW(T),S*(t-T)() H + / (X(r),(^(t-r)C)%dr 



where derivative ()' in the last term is taken over r. 
Since S* (0) = /, we obtain 



(x(t),c) H = (S(t)x ,c) H + / (S(t-r) 



A(t - a)X(a)da + A(0)X(t) 



+ / (S(t-T)*{T)dW(T),QH+ / (S(t-r)X(r),0 H dT 



X)udr 



(18) 

JO JO 
for any £ G i?*. 

Let us analyze the right hand side of Q18|). From the properties of convolution and 
(S3), 



I S(t - T)X{r)dT = - [ S(r)X(t - r)dr = - I I A{ 
Jo Jo Jo Uo 



t - s)dU(s) 



X(t - r)dr 



A(t - s)dU(s) + A(0)S{t) 



X(t - r)dr. 



We have 



and, from (S2), 



f A{0)S(t-T)X{T)dT = [ A(0)S(r)X(t-T),dr 
Jo Jo 



A(T-s)dU(s)= A(T-s)U'(s)ds= A{t - s)S(s)ds. 



Hence 



S(t-r) 



A(t - a)X{o)da 



dr= / S{t - t){A* X)(j)dr 



(S*(A*X)(T))(t) = ({S*A)(r)*X)(t) 



(S * A)(r)X(t - r)dr 



S{t)A(t - s)ds 



X{t-r)dT. 
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So, the right hand side of ifTSjl reduces to /o(5(t — r)\£ '{r)dW '(r) , if ; what means that 
(fTfijl holds. ■ 

Theorem 2 Assume that A € Lj([0,T];S ((?,#)) and i > 0, is a resolvent to the 

equation If^> £ 7V 2 (0, T; L®), then the stochastic convolution fulfills the following 
equation 

W*(t)= [ A(t - t)W* (t) dr + [ *(r)dW(r). (19) 
Jo Jo 

Proof Let us notice that the process has integrable trajectories. Then, from Dirich- 

let's formula and stochastic Fubini's theorem 



I A(t - t)W^ '(t) dr 
Jo 



(z := t — a) 



(A*S)(t — a)^{a)dW{a) 



S(t-a)^(a)dW(a) 



A(t - t) / S(t - a)^(a)dW(a) dr 
i) Jo 
t v rt 

' ^(a)dW(a) 



[ [ A(t-r)S{T-a)dT 

Jo Ua 

/*t r ft— <j ^ 
/ / Ait — o~ — z)S{z)dz 
Jo VJo 

[S{t - a) - I\${a)dW{a) 
V(a)dW(o-) . 



V{p)dW{a) 



What gives the theorem. 



4 Case A(t) = a(t)A 

In this section we consider the case when A(t) = a(t)A and, in contrary to the previous 
assumptions, A is a closed unbounded operator in H with a dense domain D(A). Then the 
equation Q has the form 

X{t)=X +[ a(t-T)AX(f)dr+ f *(r)dW(r), t€[0,T\. (20) 
Jo Jo 

For details, precise setting and several results concerning (|2()[) we refer to (Karczewska, 
2005). Because the results analogous to those from sections E] and |3] are formulated there, 
we do not repeat them here. 

We denote ^ 

W*(t):= f S(t-T)^(T)dW(r), te[0,T}; 
Jo 

now, the operators S(t), t > 0, become resolvents (see (Priiss 1993) for definition and 
properties). 

In this case, the equality (fT§|) reads 



W*(t) = [ a{t-T)AW q '(T)dT+ [ ¥(<r)dW(T). 
Jo Jo 



(21) 



Now, although the equation ()2()[) is formally simpler than lfl|). the case is interesting because 
A is unbounded. In consequence, we are not able to obtain the equality (f2*T|) directly like 
in Theorem El 
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Proposition 5 If V € N 2 {0,T;L%) and *(-,-)(J7 ) C D(A), P-a.s., then the stochastic 
convolution fulfills the equation 

(W*(t),£) H = [\a(t - T)W*(r),A*t) H + f <£, ^{r)dW(r)) H , P - a.s., 
Jo Jo 

for any t G [0, T] and £ G D(A*). 

Comment Proposition states that in case A(t) = a{t)A, a mild solution to \2ty) is a 
weak solution to \2U\) . 

This is well-known (see, e.g. (Priiss 1993)) that the scalar equation corresponding to 
(H3) is 

s(t) + n(a*s)(t) = l. (22) 

Definition 4 We say that function a G L 1 (0, T) is completely positive on [0, T] if for 

any n>0, the solution to \22\) satisfies s(t) > on [0,T]. 

Proposition 6 ( Clement and Nohel 1979) Assume that a G ^ 1 (0, T) and a is nonnegative 
and nonincreasing on [0, T]. Then a is comletely positive on [0,T]. 

Comment If a G L 1 (0,T) and is completely monotonic on [0, T], i.e. (—l) k a^ k \t) > 0, 
t G (0,T), k = 0,1, . . ., then a is comletely positive on [0, T]. 

Proposition 7 Assume that A is m-accretive operator and function a is completely pos- 
itive on [0,T]. Iffy and Afy belong to N 2 (0,T;L%) and in addition fy(-,-)(U ) C D(A), 
P-a.s., then the equality h21\) holds. 

Proof The operator A is m-accretive if an only if it generates the linear strongly con- 
tinuous semigroup of contractions, see e.g. (Zheng 2004). So, we may use results due to 
(Clement and Nohel 1979). 

Denote by A\ := j(I — J\), where J\ = (I + XA)^ 1 , A > 0, the Yosida approximation 
of the operator A. By S\(t),t > 0, we denote the resolvent operators corresponding to 
the Volterra equation (|20|) with the operator A\ instead of the operator A. The paper 
(Clement and Nohel 1979) provides the convergence lim^_^ + S\(t)x = S(t)x for t G 
[0,T] and x G D(A). 

Let us recall, that the formula ()21j) holds for any bounded operator A. Then it holds 
for operators A\, too: 

Wf(t)= f A x Wf{r)dr+ f fy{r)dW{r), 
Jo Jo 

where 

Wf{t) := / S x (t-r)fy(r)dW(r). 
Jo 

Because A x x = AJ x x for x G D(A), then 

A x Wf(t) = A\ f S x (t - r)fy(r)dW(r) = J x f S x (t - r)Afy(r)dW(r). 
Jo Jo 
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Now, from the properties of stochastic integral and by the Lebesgue dominated con- 
vergence theorem, we obtain 

lim sup E|Wf (i) - W 9 (t)\ 2 H = 
A -> 0+ te[o,T] 

and 

lim sup E\AxWx(t) - AW*(t)\% = 

what gives the required result. ■ 

Definition 5 Suppose S(t), t > 0, is a resolvent for h2U\) . S(t) is called exponentially 

bounded if there are constants M > 1 and w G R that \ \S(t)\ \ < M e wt for all t>0. 

In contrary to the case of semigroups, resolvents (if they exist) need not to be expo- 
nentially bounded even if the kernel a belongs to Existence of such resolvents is 
given, e.g. by Theorem 1.3 in (Priiss 1993). An important class of kernels providing ex- 
ponentially bounded resolvent are a(t) = t a ~ l /T{a), a G (0,2) or the class of completely 
monotonic functions. 



Proposition 8 Assume that A is the infinitesimal generator of a Co-semigroup which is 
exponentially bounded and the kernel a is a completely positive function. If \& and A^f 
belong to 7V 2 (0, T;L[>) and in addition ^f(-,-)(Uo) C D(A), P-a.s., then the equality 121]) 
holds. 



The proof bases on the following auxiliary result. 

Lemma 1 (Karczewska and Lizama 2005) Let A be the infinitesimal generator of a Co- 
semigroup T(t), t > 0, satisfying \\T(t)\\ < Me ut with M > 1, u G R and suppose 
that a e Ll QC (R + ) is a completely positive function. Then (A, a) admits an exponentially 
bounded resolvent S(t). Moreover, there exists bounded operators A n , such that (A n ,a) 
admits resolvent family S n (t) satisfying \\S n (t)\\ < Me wot , ujq G K, for all t > 0, n G N, 
and S n (t)x — > S(t)x for all x G X, t > 0. Moreover, the convergence is uniform in t on 
every compact subset o/R + . 
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